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SUMMARY 


This Report presents the results of a study of seven flows involving 
the interaction between a normal shock wave and a two-dimensional turbulent 
boundary layer. The measurements were made at free-stream Mach numbers of 
1.3, 1.4 and 1.5 and at Reynolds numbers based on an effective streamwise 
run of 10 ~« 10° to 30 x 10°, The results were obtained from comprehensive 


traverses with both pitot and static probes. 


Standard boundary-layer integral parameters based on wall and measured 
static pressures are presented, together with velocity profiles and the Mach 


number distribution over the interaction region. 


An investigation has been made of the 'law of the wall' and the ‘law 


of the wake' under the influence of strong normal pressure gradients. 
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AANZY 


1 INTRODUCTION 


A series of experiments!” has been undertaken in the RAE Bedford 3ft x 3ft wind 
tunnel to investigate the interaction of a normal shock wave and a turbulent boundary 
layer at nominal upstream Mach numbers of 1.3, 1.4 and 1.5 over a Reynolds number range 
based on the undisturbed boundary-layer momentum thickness at the start of the inter- 


action 6£°12.% 10° to 34° 10". 


This Report deals with the measurements made at a series of stations upstream and 
downstream of the interaction using conventional pitot and static probes traversed normal 


to the flat wall. 


The series of experiments was planned because of the lack of knowledge of the inter- 
action between normal shock waves and turbulent boundary layers at high Reynolds numbers 
and the need to predict the development of the boundary layer through and beyond the 


interaction region as, for example, in the flow over supercritical aerofoils. 


A number of investigations have been reported over the past two decades and broadly 


speaking three techniques have been used to produce the interaction. 


The first and most commonly used technique (Fig 1a) in supersonic tunnels is to 
position a shock generator with downstream choking flap over a flat plate so that the 
steady normal shock wave formed by the shock generator interacts with a turbulent 
boundary layer grown from the leading edge of the flat plate. This is the technique 


adopted by Seddon”, Vidal, et at® and Kooi’: 


The second technique’? 


is used in transonic wind tunnels and involves mounting a 
two-dimensional bump on the wind tunnel wall so that the shock-wave boundary-layer inter- 
action approximates to that on an aerofoil (Fig 1b). More recently, small supercritical 
aerofoils have been tested! ° and in fact these two experiments have been combined by 
Burdges!! who let the aerofoil into the tunnel floor and bled away the floor boundary 
layer under its leading edge. 


Gada! ? and more recently Mateer, Brosch and Veigas!? used a third technique which 


was to hold the normal shock wave steady in a supersonic tube by the adjustment of a 
conical choke downstream of the working section as in Fig Ic. The turbulent boundary 


layer under consideration was allowed to grow naturally along the tube wall. 


For the present experiment the facility which was available was the RAE 3ft « 3ft 
transonic~supersonic wind tunnel, The arrangements of Fig la or Ib could have been 
employed, but it was argued that if the tunnel floor was used, then the measuring tech- 
niques would be simpler, a higher Reynolds number could be obtained, and the inter- 
ference effects of the wall boundary layers should be less than in the technique of 
Fig la where the test boundary layer is of smaller thickness than the interfering wall 
boundary layers, The simple technique was therefore used of holding a shock wave across 
the test section, set up to run at supersonic speed, by means of an adjustable sonic 
throat which was situated far downstream of the test section. The technique is in fact 


similar to that of Fig Ic. 


The 3tt + 3ft wind tunnel has a good Reynolds number capability (up to 12 ~ 10°/m 
tor continuous running) is easily accessible for non-intrusive measurements and, above 
all, is easy to modify i-, the working section region having removable supersonic liners 
and downstream wooden fairings. A further advantage is that the pressure distribution in 
the working section is not only similar to that over a supercritical aerofoil but also 
the dimensions and Reynolds numbers approximate to full-scale conditions. A ninth scale 


(4 in « 4 in) model of the tunnel is also available and this proved invaluable for the 


development of the experiment. 


the most time-consuming aspect of the development has been providing a suitable 
downstream sonic throat to control the shock-wave position and to keep it acceptably 
steady. However it was also necessary to design and manufacture a new raised false tunnel 
floor (to house submerged boundary-layer traverse mechanisms and to raise the boundary- 
layer interaction region into view above the bottom of the schlieren windows) which would 
match the three alternative half-nozzle blocks used to generate the required free-stream 


Mach numbers and which form the upper wall. 


During the period of time taken to manufacture the new floor and the necessary 
associated equipment, two interim experiments were made, the first” using a conventional 
boundary-layer pitot rake and the second! using non~intrusive laser-Doppler measurements. 
The first experiment differed in another important aspect from the present one in that 
the pitot rake was attached to the tunnel floor in a single, fixed, position and the 


shock wave moved fore-and-aft of the rake. 


For the present experiment the location of the shock wave was fixed for each condi- 
tion tested and measurements were made with probes which could be traversed both normal 
to the tunnel floor and streamwise. Measurements of the static pressure distribution 


along the floor were also made. 


7 DESCRIPTION OF THE EXPERIMENT 


deot Mechanical arrangement 


The general arrangement of the experiment is shown in Fig 2. The unmodified parts 


of the 3ft = 3ft wind tunnel are shown shaded, 


In its standard form, the 3ft <x 3ft wind tunnel has a working section with a fixed 
lower Liner and interchangeable upper liners which are matched to the lower liner to 
give 0.1 incremental steps in Mach number between 1.3 and 2,0, One of the major modifi- 
cations made to the tunnel was the provision of a new raised bottom liner together with a 
false floor (Fig 3). This false floor contained two pitot traverse mechanisms so that a 
region close to the floor centre line could be investigated from approximately 0.8 m 
ahead of the schlieren window centre line to 3m downstream using a total horizontal tra- 
verse movement of little more than 2 m. Because of the interest in the free-stream Mach 
number range of 1.3 to 1.5, the new bottom liner was designed to match the M = 1.4 upper 
liner and the consequent slight mismatch with the M = 1.3 and M = 1.5 upper liners had 


to be accepted. Although raising the floor by 152 mm had the advantage of bringing the 


a 


shock-wave interaction region into full view through the schlieren windows, it reduced 


boll 


the height of the working section so that the tunnel was no longer square in cross 


section. There was however no particular disadvantage in this because the asymmetrical 
nozzle which generated the supersonic flow inevitably resulted in boundary layers which 


were not the same on al] four walls. 


The other major item of redesign as described in Ref 2 involved the careful fairing 
of the existing 'spoiler door’ arrangement* in the diffuser section to give an adjustable 
second throat. By adjusting the throat, the normal shock wave could be placed in the 
desired location wi-hin the view through the tunnel windows. Any subsequent small pro- 
gressive movement of the shock wave could then be corrected by adjusting the tunnel 
compressor speed. Initially there was a large random movement of the normal shock wave 
but this was reduced to about 20 mm by the careful fairing of the shape of the second 
throat. The mean position of the normal shock wave was checked both visually and by 
means of a differential pressure transducer connected between a wall static pressure hole 


under the shock wave and a reference tapping near the working section sonic throat, 


The region between the working section and diffuser was faired to give smooth side- 
wall and roof contours. The geometry of the fairing was such as to give an expanding 
passage, the cross sectional area at the diffuser entry increasing to 1.067 times the 
area of the working section over a distance of 3660 mm. In the region of the interaction, 
the standard taper of 0.004 m/m was maintained on the top liner. This is normally com- 
bined with a similar taper on the lower liner to allow for boundary-layer growth on all 
four tunnel walls. However, for this experiment, the taper on the false floor was 


approximately 0.003 m/m. 


The traverse mechanism together with a twin pitot probe is shown in Fig 4. Two 
identical traverse mechanisms were supported 1800 mm apart from a slide which was 
approximately 6 m long and sufficiently flexible to submerge, at its upstream end, below 
the liner surface approximately 1 m ahead of the schlieren window centre line, and at its 
downstream end to retract into a covered channel between the spoiler doors in the diffuser 
section. Any leaks around the slide were sealed by two pvc tubes which inflated when the 
slide was stationary as shown in Fig 3. A chain drive was used to move the slide in the 
streamwise direction. A total streamwise distance of 3800 mm could be covered using the 
tandem boundary layer traverses. The slide was positioned 100 mm to the port of the row 
of static pressure holes on the tunnel centre line and the probes were cranked so that 
vertical traverses could be made half way between the slide and the tunnel centre line. 
In this way it was hoped to reduce the various interferences, namely, of the probe on the 
wall static pressure holes, of the slight irregularicy of the slide on the probe measure- 


ments, and of the probe stem on the probe measurements. 
Repeatability of the probe position was +0.02 mm vertically and +3 mm horizontally. 


Details of the static probes can be seen in Fig 5. The pitot probes were designed 
to be fairly short in order to avoid vibrations in the expected highly turbulent flow. 
The same overall dimensions were then retained for the twin static probes. With hind- 
sight, the twin static probes might have been lengthened to diminish the interference 


EE I 


* Used in normal testing for controlling the tunnel shock wave during starting or 
stopping the tunnel, 


caused by the supporting structure because the vertical positional accuracy is not quite 
as important for these probes as for the pitot probes. However this was not done and 
large corrections had to be applied to the static pressure measurements especially at 


transonic Mach numbers, These corrections are described in the Appendix. 


All the probes were electrically insulated from their vertical traverse mechanisms 
so that a touch indicator could be used to set the pitot datum at the surface of the 
talse floor. The streamwise datum was set by aligning the tip of a pitot probe or the 


holes in a static probe with a line marked on the surface of che false floor. 


Pressure measurements were made with differential transducers (Druck) of range 
“69 kN/me in two D-type scanivalves placed outside the tunnel shell. The transducers were 
calibrated against a Texas Instrument 0-203 kN /m” absolute quartz Bourdon-tube pressure- 
controller. Corrections were made to the primary slopes and zeros of these calibrations 
for each data point, by comparing the transducer outputs with reference pressures applied 


to the first and last pairs of scanivalve ports. 
2.2 Experimental measurements 
2.2.1 Pitot measurements 


Measurements were made at three nominal Mac numbers, 1.3, 1[.4 and 1[.5 at Reynolds 
numbers of approximately 10 x 10° /m (the maximum continuous running value permitted for 
the experiment) and 3.5 ~x 10°/m. An intermediate Reynolds number was also included for 


M= 1.5. Details of the test conditions are given in Table }. 


Traverses normal to the floor were made using the twin-pitot probe over a large 
number of stations f om 0.7 m ahead of the normal shock wave to approximately 3 m down- 
stream of the normal shock wave, In most cases, a vertical distance of 200 mm was 
covered, The forward pitot was removed while measurements were being mace with the rear 


one, and a fiush plug was used to fill the hole left in the slide. 


-2 Floor static pressure measurements 


Measurements of the static pressure distribution on the floor were made with only 
the rear pitot in place, set at its furthest downstream position. These measurements 
were used to calculate the wall reference static pressure for each pitot or static 


pressure probe traverse, 


2.2.3 Reversed pitot measurements 


Measurements were also made with a single reversed pitot probe (rather similar in 
shape to static probe B (Fig 5) but reversed) for conditions of separated flow which 
occurred for M = 1.5. Separation was indicated by the forward-pointing pitot recording 
a lower pressure than the corresponding wall static pressure hole. The measurements 
were made with the tip of the reversed pitot in the same streamwise position as the 
normal pitot tube. Although the flow was probably disturbed by either probe, a simple 


correction was made to both sets of results as described under section 3.1, 


2.2.4 Static probe measurements 


Before the static pressure measurements in the boundary layer were begun, the 


g 


v 


660 


static probes were calibrated in the centre of the slotted transonic working section of 


the tunnel over a Mach number range of 0.4 to 1,2 and similar free~stream Reynolds 

numbers to those encountered in the main experiment. Each probe was held at the end ot 
its stem furthest away from the probe tip by a specially manufactured sting. Static probe 
B was also calibrated while supported at the near end of its stem to check for support 
interference. This appeared to be negligible. The probe measurements were compared with 
wall static pressures obtained from a tunnel calibration referred to a 3° conical static 
head probe. This probe has been described in Ref 14 and within the limits of the required 


accuracy was assumed to have no static pressure errors. 


The real problem which arose when using static probes in an environment which 
invoked the problems ». transonic flow combined with pressure gradients, wall interference 
etfects and the presence of boundary layers, lay in obtaining sufficient calibration data 
in a tunnel which could nut be used petween M = 0.88 and M = 1.3 while modified to take 


the false floor and traverse mechanisms. 


However some information was forthcoming from twin static probe calibrations made 
through the floor boundary layer between M = 0.66 and 0.88 at similar Reynolds numbers to 
the main experiment. Additional information was obtained during the actual shock-wave 
boundary-layer interaction experiment from the static and pitot probe traverses ahead of 
the interaction, while further data were available from traverses | to 3 m downstream of 
the normal shock wave at M = 1.3. In all these cases it was assumed that there was no 
static pressure variation normal to the wall through each vertical boundary-layer tra- 
verse, and that the pressure was equal to the estimated wall value. Because the free- 
stream calibration showed no apparent Reynolds number effects, data at the highest 


Reynolds number from the experiment were used. 


Time limitations meant that nearly all the boundary layer static traverses were 
made with the twin-static probe rather than the single probes and the static investiga- 
tion was omitted for M = 1.4 at a Reynolds number of 3.5 x 10° /m. Single static probes 
were in fact used only as a check on the validity of the results obtained from the twin 
probe. Traverses were made in rather different streamwise positions from those made with 
the pitot probes and each traverse contained about half the measurements made using the 


pitot. 


2.2.5 Shock-wave position 


A careful check was made of the shock-wave position using the output from the 
shock-position transducer and readings were only taken while the shock wave was within 
prescribed limits. Subsequently, during the analysis of the results, the average posi- 
tion of the normal part of the shock wave was ascertained from the large number of 
schlieren photographs taken for each set of tunnel conditions. This was necessary 
because there was an interval of a year between making the static pressure measurements 
and the pitot measurements. It was difficult to reset the shock wave in the same posi- 


tion, and so corrections had to be applied to make both sets of results compatible. 
3 REDUCTION OF DATA 


The data have been reduced as follows to produce the boundary layer profiles of 


Table 2 and the integral parameters of Table 3. 


3.1 Boundary~layer profiles 


These have been tabulated against y (the equivalent height of the pitot or static 


tube above the wall) where for pitot measurements 


h+0.15d (1) 


<< 
W 


and for static measurements 


where h = height of probe centre line above floor 
and d 


diameter of probe. 
The profiles have been derived 


(a) assuming constant wall static pressure applied throughout each traverse, 


(b) using measured static pressures after correction (see Appendix). 


Standard formulae were used to calculate Mach numbers, velocities and densities 
from the pitot and static results. Total temperature was assumed constant across the 
boundary layer except at the wall where a recovery factor of 0.89 was used. Wall static 
pressures were calculated for each traverse position by linear interpolation of the floor 
static measurements made with the front pitot removed and the rear pitot as far down- 


stream as possible. 


The flow was assumed reversed when the pitot pressure was lower than the wall 
static pressure. In such cases measurements were also made with the pitct tube reversed 
in direction. The Mach number was taken as the mean derived from the two sets of measure- 
ments. When the pitot tube was effectively facing downstream, it was assumed that the 
pressure recorded was that for a base with a pressure coefficient of -0.6. The wall 


static pressure was used for both sets of calculations in the reversed flow region. 


Unit Reynolds numbers were calculated using the following formula based on 


Sutherland's law for viscosity 


MePs 3 
Re/m = = (T, + 110.4) <x 47.91 «x 10 (3) 


a 


where subscript 6 denotes conditions at the edge of the boundary layer 
p is the static pressure in N/m? 


and T is the temperature in K. 


A further parameter, namely P;/Ptg has been presented where P; is an estimated 
static pressure in the ‘equivalent inviscid flow' (see Ref 15). The static pressure Ps 
is that in the equivalent inviscid flow which is defined as the flow external to the 
shear layers continued as a: inviscid flow to the wall bounding the real flow with the 
growth of the viscous layer represented by transpiration at the boundary. P; has been 


non-dimensionalised by dividing by the tunnel total pressure Pty - 


HAOTS 


An estimate of Bil P eg has been attempted in the region ahead of the normal shock- 
wave, Here, the flow at the edge of the boundary layer is supersonic and influence lines 
of constant pipe, following Mach lines may be constructed back to the wall from the 
edge of the viscous region. An estimate of P,/Pt, as a function of y for each tra~ 
verse can then be made from the carpet of influence lines which are assumed straight and 
inclined at the Mach angle yu plus the flow inclination at 59 999 (the distance from 


the wall where U/U, = 0.999) viz 


Vv aa* d are 
where (%) = “ax 7 (So.999 7 8) ax (inc ,t,)) oa 


and Vis the vertical component of velocity 
U is the horizontal component of velocity 
6* is the displacement thickness 


o is the density. 


It will be noted (see section 4.4) that there is a supersonic region at the edge of 
the boundary layer behind the normal shock wave for M = 1.5, However estimates of 
P,/Pe,, have not been included for this region because they become discontinuous in the 


region of the shock wave, 


3.2 Integral parameters 


Boundary-layer thickness parameters were obtained by trapezoidal integration. 
Velocity profiles were faired between the wall and the Pitot position corresponding to 
i eer lo , 
0.8 mm from the wall by applying East's prediction of the law of the wall In compress~- 


ible boundary layers as a 20 point curve, 


East's prediction is based on an incompressible law of the wall combined with a 
Be du : , 17 ae 
compressibility factor. For this experiment Cole's law of the wall has been used with 


. 2 18 
the constants recommended at the Stanford conference name ly 


yu 
Ui T ‘ 
T = 5.62 log ye + 5.0 , (5) 


v 


where U_ is the friction velocity. 


The form of East's prediction which is faired into Cole's law of the wall is 


U \ : 2K -y*t/D : 
te viol Ff2x n(y* Re ') + o(: op tt )] . (4) 
T 


where the compressibility factor 


7 FOE) ot 
F M. 7 Vie x 0.021033 , 
2 
D = 8.73(1 + 45F°) , 
Ve 
y* = Vines 
Ww 
K = 0.41 , 


r being the recovery factor, and suffix w denoting wall values. 


A value of U, was obtained at y * 0.8 mm using equation (6) and values of 


U/U, were then calculated for 20 equispaced values of y between O and 0.8 m, 
A linear variation was assumed between y = 0 and 0.8 mm for reversed flow. 


The following definitions of the integral quantities have been used as first pro- 


eee | é a dors Fi F 
posed by Myring 9 and later incorporated in East's modified momentum integral eyuacivn’>: 


They are (where suffices i and w denote ecuivalent inviscid flow quantities and 


wall values respectively) 


6 
en ea [ow, -pU)dy , (7) 
0 


loc) 
i) 


QO. . 
iw lw 


§ 
es 2 2 
§ = ef f{(o,07 oU ) B, tou; ouybay ; (8) 
Piw iw 0 


shape parameter 7 (9) 


x= 
u 
DI fot 


So. U, 
Piw iw 0 


6 
shape parameter H = —— Jou, -U)dy . (1.0) 


The integral parameters were calculated for three different distributions of static 


pressure across the boundary layers as described in the following paragraphs. 


3.2.1 Static pressure constant and equal to the measured wall pressure 


Equations (7), (8), (9) and (10) reduce to the familiar standard integrals when 


static pressure is constant across the boundary layer giving 


ey, 


leva) 


F : cU ‘ 
displacement thickness se = [() 7 so ) dy , tla 
. ATA 
6 ee 
\ 
$6 
; 7 oU U - 
momentum thickness ai [ ai (: = ¢) dy , pat 
ae Tes § 
0 
o* 
shape factor H =>, (13% 


- ] fe i. 
shape factor H => [ = (: = =) dy. C14) 


& 9 
: ! a U 7 é 
The energy thickness was also evaluated: é = —— | I -(— dvr. (15) 


3.2.2 Static pressure variation as measured 
Equations (7), (8), (9) and (10) can be used in an ‘intermediate form' as in Cook. 
aie oe : , : 
McDonald and Firmin OL In this form the boundary layer defect thicknesses are refurred 
to a fictitious potential flow having the same static pressure distribution as that for 
the actual viscous flow and a total pressure equal to that at the edge of the boundary 


layer (pty) 


Thus for example 


(lAD 


where Pr 


The revised integral quantities are as follows: 


6 
displacement thickness 6* = [G U_ - pU)dy , (18) 
fy) 


6 
‘ 2 2 
momentum thickness Os 2= i {(° UT = pU } - U(e_Uo = cl bay (19) 
PP pp p cU)pdy , 


, | 
shape factor H = zr: 
7 | 
shape factor H = eee be (U. -U)dy . (20) 
P 6poU J op op 
Woe { 


3.2.3 Equivalent inviscid static pressure together with measured static pressure 


Aq 


» =, H and H have already been defined by equations (7), (8), (9) and (10). 


Equivalent inviscid static pressures were estimated as described in section 3.1 for 
traverse positions ahead of the normal shock. Equivalent inviscid Mach numbers were 
obtained from these pressures, and the total pressure calculated at the edge of the 


boundary layer (equation (17)) so that 


: (21) 


cy and Us were obtained from M; using a temperature recovery factor of unity 
while the values of ge and U were identical to those in section 3.2.2 (being calculated 


from the measured pitot and static pressures). 


Linear interpolation was used to obtain Oo; and U; for the 20 points between the 


wall and y 20.8 mm. 
3.3 Skin friction 
Three estimates of skin friction are given. 


The first estimate ES has been deduced by considering the pitot tube when in con- 
tact with the wall to be a Preston tube and applying Patel's*! calibration as formulated 
by Head and Vasanta Ram” and transformed for compressible flow by the method of Fenter 
and Seainack =. 

The second and third estimates emerge directly from ul which was calculated by 
fitting East's law of the wall (equation (6)) to the pitot reading for y * 0,8 mm. 

The second estimate was obtained assuming constant wall static pressure to apply across 
the boundary layer and the third was obtained using measured static pressures. Because 
at y * 0.8 mm wall static pressure applied for both cases, the differences are entirely 


due to the boundary layer edge conditions obtained from the different static pressures. 


It will be noted that several values of Cg are missing from Table 3 and these 


omissions occur where reversed flow is indicated at M = 1.5. 


4 EXPERIMENTAL RESULTS 


4.1 The general characteristics of the flow 
4.1.1 Shock-wave patterns 


Fig 6a shows the vertical cross sections of the main interactions traced from 
schlieren photographs. Typical photographs have been reproduced in Fig 6b. For complete- 


ness, the calculated boundary layer thickness (6 ) and the start of the measured 


0.995 
reversed flow regions at M = 1.5 have been added to the diagrams. The progression from 
the pattern at M = 1.3 of a single normal shock wave distorted by the influence of com 


pression curves emanating from the thickening of the upstream boundary layer, to the 
established triple shock system with the shear layer emanating from the point of bifurca- 
tion at M = 1.5, is well illustrated. The increase in size of the shock system with 


reduction in Reynolds number can also be seen. 


The main interaction region has been described in detail by East? and the present 
results at high Reynolds numbers agree with his at M = 1.3 and !.4. However there is a 
discrepancy between the sets of results at M = 1.5. The shock bifurcation point for the 
present experiments is at 180 mm from the floor compared with a value, presumed to be 


erroneous, of 215 mm given by East. 
4.4.2 Flow visualisation 


After the completion of the traverse measurements the interaction region was 


examined by oil flow*, 


Photographs of the surface oil flow were taken for all three Mach numbers at a 
Reynolds number of 10 ~x 10°/m. These photographs, reproduced in Fig 7, were taken after 
shutting down the tunnel and removing the roof and therefore suffer from slight blurring. 
The viewpoint looks downstream, The photographs reveal mild three-dimensional effects 
for M = 1.3 and M = 1,4 but a much more complex pattern for M = 1.5 when the flow is 
separated. The photograph for M = 1.5 is shown in more detail in Fig 8 together with an 


attempted interpretation of the wall streamlines. 


The nature of the flow in this case appears to be such that there is flow along the 
arch of a vortex connecting a position on the port side, denoted by B in Fig 8, tw a 
position denoted by A on the starboard side. Thus in a sense A_ is an attachment node 
and B a separation node. It should be remarked that whilst the streamwise ex-ent of 


this region is roughly 300 mm the probe measurements suggest its depth to be only 7 mm. 


Near the tunnel centre line between the two vortex patterns, the flow is tolerably 
two-dimensional except in the immediate regions of the saddle points at the start and 


end of the separation region. 


Presumably because of the non-uniform boundary-layer thickness on the tunnel side 
wall arising from the flow field associated with the assymetrical nozzle, the separation 


line on the sidewall is swept and a single vortex node only of separation type is formed. 


* The oil-flow mixture was an amalgam of the following in the ratio of 4 cc Vitrea Fl, 
to 2 cc Limea 931, to 3 cc Ti0,, to 2 drops of oleic acid. 


14 


The vortex arising from this node is presumed to create the strong convergence towards 
the tunnel centre line occurring in the corner regions downstream of the interaction. 
The lines D in Fig 8 have the appearance of separation lines but may in fact indicate 


only a locally strong convergence of the boundary layer. 
4.1.3 Momentum balance 


Momentum balance calculations have been made using a rearranged form of East's 


ee , noel 
sodified momentum integral equation 


a a dX 
iw iw 
= uM? dv, av 
iS us d Ko. v2 3 H = J] + —1W__iw lw : (22) 
2 dX lw iw = 2 dX 
p. U. H U. 
iw lw iw 


The last term has been neglected in the analysis because it is smaJl and in any 
case cannot be evaluated with any reliability. The remaining terms in the right hand 


side of the equation are: 

@ the normal pressure gradient effect arising from the wall and stream curvature; 
Q the approximation to the direct effect of wall curvature on the flow momentum; 
© the approximation to the normal stress terms, 


The integral parameters, H, i, 6 and 6§ are defined in equations (7) to (10) 


iw 
horizontal component of velocity and Mach number at the wall. 


while ony U, and Mi are respectively the equivalent inviscid values of density, 


The surface curvature « is zero for this example leaving 


48 Panta ee hy 
a (H+ 2-H) ge 2 
lw ep. U. 
1w lw 

1 d Pa ay 1 ‘| 2 <-fH - 1 

= = 10. Ut 8) + — (Ko. UT i( = )f (23) 
20. v2 dx } iw iw ax? et y? dX iw iw i 
lw lw lw 1W 


3 . 2 . ; 
After multiplying throughout by Oey and rearranging, the equation becomes 


ce 


6bo0¢ 


¢ 


(24) 


2 
Hat) 


f = 0.072( = » and is the contribution due to the normal stress, 
H 
eae 2 aes ~_¢ 
and f, Sy aa = (H + 6°) = a (1 + H)(6 + 6) and is the contribution due to the 
dx” H ax2 stream curvature 


If suffix 0 represents the conditions at the start of the interaction, then 


dividing through by Piwqliwg’ ig and integrating, results in the form of the momentum 


integral equation used in the momentum balance calculations shown in Fig 9: 


© (tx) © 


is iw = es 
Gee 8 ae oi OF tee ee agri IN fp. Ur 8 
iw iw -] [4 ee ax ele [ iw 4/—iw a iw iw ¢ 
2 = 2 2 2 mie i 262 0) 
Piwqliwg®o ky Piwqliwg®o 1 Piwg®o *¥O Piwqliwg®o 
sragesereey C25) 
P iWin! 
In Figs 9 and 11 |——-—-;-— - 1 is called the left hand side while the remainder 
PiwQiwg’g 


of the expression is called the right hand side. The numbered terms are identified later 


(in Fig 12) as the contributions due to 


@ skin friction 


® pressure gradient. 


Both sides of equation (25) were calculated using the alternative assumptions of 


section 3.2 which were (for traverses across the boundary layer): 


(a) static pressure constant and equal to the measured wall pressure; 
(b) static pressure as measured; 
(c) equivalent inviscid static pressure together with measured static pressure. 


The equivalent inviscid static pressure was calculated for traverses ahead 


of the shock wave only. 


Fig 9 shows the comparison between the left hand side and the right hand side of 
equation (25) plotted against the streamwise position. The calculations were made using 
assumption (c). It will be seen that the momentum balance is good over the main inter- 
action region and only becomes significantly in error at approximately 800 mm downstream 
of the normal shock wave. Downstream of this point, the left hand side actually reduces 
fo: Mach numbers ahead of the interaction of 1.4 and 1.5. This indicates a flow diver- 


gence and it is interesting to calculate its magnitude. 


In terms of equation (22) Green et a" showed that a flow divergence could be 
accounted for by subtracting the local rate of divergence 6(d¢/dz) from the right hand 
side. This quantity represents the rate of the deviation relative to the nominal flow 


direction with respect to distance Z perpendicular to the flow but parallel to the surface. 


In equation (25) this divergence takes the form 


xX u2 = 
Piwiw do ax 
. yt 9. a2 
xX e iw’ 1Wo fe) 


and this can be added to the left hand side. It is therefore easy to make an estimate 


of do/dz 


This estimate is shown in Fig 10a-c in degrees per metre. For points further down~- 
stream of the normal shock wave than 300 mm the figures show a flow divergence of 7° /m 
at M = 1.3 rising to 15°/m at M = 1.5 (equivalent to 7° and 15° over the tunnel width). 
The wild fluctuations within the main interaction region are of course due to the 


rapidly changing conditions and the close streamwise spacing of the traverse positions. 


Assumptions (a), (b) and (c) were used to provide comparisons of both sides of the 
momentum integral equation as shown in Fig 11. There appears little to choose between 
the momentum balance using any of these assumptions, However it is possible to see some 
improvement using assumption (b) rather than (a). If the region ahead of the shock wave 


is examined there is again some improvement if calculation assumption (c) is used. 


Fig 12 shows the individual contributions to the right hand side of equation (25). 


4.2 Edge Mach numbers and boundary layer parameters 


Fig 13a-c shows the development of the Mach number at the edge of the boundary 
layer, M, , the boundary-layer displacement thickness, ¢* , the shape factor H, 


defined by equation (13), the boundary layer thickness, the skin friction 


_ °0.995 ? 
coefficient ae (section 3,3), and the shape factor H defined by equation (14), all 
plotted against the longitudinal distance X from the normal part of the shock wave. 
The results in Fig !3 were calculated assuming constant static pressure, Pp, » across the 
boundary layer normal to the tunnel wall. Further graphs of Mas §*, H and H 
(Fig 14a-c) show the effect of using measured static pressures to calculate these para- 
meters. The longitudinal region covered by these figures is from 500 mm ahead to 500 mm 
downstream of the normal shock wave. Outside this region it was assumed that the static 


pressures were constant normal to the wall and equal to the static pressures used in the 


first set of calculations (see Fig 15a&b and the Appendix). 


Overall there appears to be little effect of variation of Reynolds number on the 


shape factors and the normalised variation of 6* and A separation bubble 


°9.995 7 
occurs under the interaction region at M = 1,5 and its extent is indicated by the zero 
values of ff . The bubble size decreases slightly with increase in Reynolds number. 


The shape of the bubble can, however, be seen in the plots of Mach number distribution, 


AHOCH 


Fig l6a-f, and will be described in section 4.4, 


The results agree broadly with those obtained in Ref 2 where a fixed pitot rake was 
used on the standard tunnel floor. However discrepancies begin to creep in | m downstream 
of the normal shock wave where in Ref 2 the Mach number at the start of the interaction 
was atfected by the proximity of the tunnel throat. The results from Ref 2 at M = 1.5 


were much more sensitive to Reynolds number in the interaction region and this was 


probably due to pitot-rake interference. One other point should be noted (and it can be 
seen in nearly all the following figures) and that is that the rehabilitation process 
behind the normal shock wave has not been completed after 3m. This distance is equiva- 
lent to at least 400 times the undisturbed displacement thickness ahead of the inter- 
action, 1t can be seen in Fig 13a-c that the edge Mach numbers, shape factors and 
boundary layer thicknesses have not yet reached their asymptotic values at the furthest 


downstream traverse positions, 


4.3 Measured static pressures 


Static pressure distributions normal to the wall are plotted on scale diagrams of 
the main interaction region in Fig 15a&b as the difference from the corresponding wall 
static pressure non-dimensionalised by dividing by the undisturbed stagnation pressure. 


Zero values are located at the appropriate streamwise station. 


As described in the Appendix, the static pressure results have usually been made to 
merge into the appropriate wall values at floor level, however wall static values have 
been assumed to apply throughout the separated region in the absence of reliable local 
values. Static pressures recorded very close to the shock wave may have been affected by 


its movement and so local pressure gradients may have been reduced. 


The diagrams show the expected static pressure variation throughout the main inter- 
action region. The static pressures are fairly constant across the inner part of the 
shear flow and then decrease outwards across the compression region ahead of the normal 
shock wave to free-stream values. Behind the main shock wave they increase from one 
fairly constant level near the wall to a second level behind the normal part of the 


shock wave. 
4.4 Mach number distribution 


Mach number contours are plotted in Fig l6a-f over the main interaction region and 
over the whole of the region investigated downstream of the normal shock wave (note the 
change of streamwise scale). Measured pitot and the corrected static pressures were used 
to calculate Mach numbers and the figures result from linear interpolations of Mach 
number profiles in the vertical direction and from horizontal plots of Mach numbers at 


constant distances from the wall, y. 


The results are therefore dependent on the accuracy of the corrected static 
pressure measurements and the interpolation procedures which are affected most critically 
by the smearing effect of the slightly unsteady shock wave. However the results in the 
interaction region are broadly confirmed by the laser-Doppler anemometer measurements of 


East? and the results of Kooi! at M= 1.4, A particular feature is the development of 


the supersonic tongue beneath the downstream shock wave for M = 1.5. 


If East's and the present results at a Reynolds number of 10 ~« 10°/m are compared 
in slightly more detail, the Mach number contours are in good agreement ahead of the 
normal shock wave if due allowance is made for the smearing effect of the unsteady shock 
wave (see above). Very similar supersonic tongues may be seen behind the shock wave and { 
they continue downstream (in both experiments) for !00 mm near the edge of the boundary 
layer. East's results at M = 1.5 do not cover the full extent of the supersonic tongue 
but show that it is more extensive than at M = 1.4, In fact the present experiment con~ 
firms that it persists near the edge of the boundary layer for 400 mm downstream of the 


shock wave while at a Reynolds number of 3.5 ~< 10°/m it persists for 600 mm. 


It is difficult to be specific about the shape of the supersonic tongue immediately 


gk : of 
downstream of the trailing shock wave, but the present results agree well with Kooi who 
also used pitot and static pressure measurements, The rear shock wave appears to termin- 
ate at the edge of the boundary layer so that supersonic flow appears to exist above the 


boundary layer behind a normal shock wave, 


There is evidence to suggest that the thin shear layer formed downstream of the 
shock bifurcation point is still intact 3 m behind the normal shock wave, and even at 
M = 1.3 there is no sign that the inviscid flow has recovered sufficiently to provide a 


uniform Mach number contour across the tunnel. 


The other point of interest is the separation region at M = 1.5, There appears to 


be a small decrease in the overall dimensions of the separation bubble with increase in 


Reynolds number. In fact the length of the separation bubble is reduced from 300 mm at 

Re/m = 3.5 x 10° to 200 mm at Re/m = 10 ~« 10° while the height is reduced slightly from 
just over 9 mm at the lowest Reynolds number to 7 mm at the highest Reynolds number. This 
result is in disagreement with Sawyer et aie where a much greater reduction in bubble size 


was indicated although as stated in section 4.2 the conditions were somewhat different. 


4.5 Boundary~layer velocity profiles 
4.5.1 General 


Boundary-layer profiles are shown in Fig 17a-f plotted on a scale diagram of the 
main interaction region together with an extended diagram on a compressed scale showing 


the whole region of measurement behind the normal shock wave. 


The figures indicate the general characteristics of the boundary-layer development. 
They show the progression from the undisturbed state through the leading compression 
which influences the shape of the outer part of the boundary layer, through the subse- 
quent separation or near-separation region, to the final recovery towards a zero-pressure- 
gradient form for the profiles. It can be seen even on the small scale diagrams that 
the profiles are still disturbed after a distance of nearly 3 m downstream of the normal 


shock wave, or approximately 400 times the undisturbed displacement thickness. 


4.5.2 Logarithmic velocity profiles 


Representative velocity profiles are plotted in the logarithmic form adopted by 


; 2 2 . i i . . 
Winter and Gaudet : as (urut) against log(yUt/v5) where ur 1s the equivalent 
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incompressible friction velocity. Values of ur were derived from the skin friction 
coefficient CE. using the formula 

h 
i “fp } 2 2)" (26) 
vu = us 51 ( + 0. M. . 


Ce was obtained by treating the pitot tube when in contact with the wall as a 


Preston tube and using the calculation method described in section 3.3, 


The results are plotted separately for three regions. Fig 18 shows the undisturbed 
velocity profiles just ahead of the interaction region. These profiles were used to 
define the undisturbed conditions of Table 1. The profiles far downstream of the inter- 
action region are shown in Fig 19. Fig 20 shows typical profiles just downstream of the 


main interaction region. 
4.5.3 The law of the wall 


Figs 18 to 20 show typical ‘law of the wall' fits for small values of y (the dis- 
tance from the wall). The law of the wall is defined in the form adopted by Winter and 


Gaudet”? as 


(27) 
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which may be expressed as 


(28) 


fa} 
a wis 
W 
> 
~— 
ie] 
a 
= 
MS 
ee a 
m 
“ 
+ 
w 


The results have been compared with Winter and Gaudet's incompressible law of the 


wall 


= 6.05 108 (=) + 4.05, (29) 
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Coles’ incompressible law of the wall is usually applied to compressible flow by 
making use of the Van Driest~° transformation. It was therefore necessary to confirm 
that the transformations of Van Driest and Winter and Gaudet gave similar results when 
applied to the compressible profiles of the present experiment. Two profiles were 
checked and are plotted using both transformations in Figs 18 and 19, The first profile 


had a free-stream Mach number of 1.54 and was in the undisturbed region ahead of the 
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interaction region, while the second profile had a free~stream Mach number of 0.83 and was 
from a pitot traverse made 2759 mm downstream of the normal shock wave. It can be seen 
that both transformations gave identical values of A and B although the wake compo~- 
nents were very different at the free-stream Mach number of 1.54 (Fig 18). Because both 
Lransformations gave identical values of A and B Coles‘ incompressible law of the wall 


has been presented in the form adopted by Winter and Gauder, 


It would appear from Fig 18 that within the data scatter there is little to choose 


between the forms of the law of the wall due to Coles and due to Winter and Gaudet, 


The average fit made to all the profiles ahead of the interaction region is 


at 
U y tT 
— = 6.16 log — + 3.6 (30) 
vr Ys 


which is somewhat closer to the latter. 


However behind the normal shock wave (Figs 19 and 20) the average fit becomes 


i 


yu 
Oat RTE Sterns oy thee (G31) 
ur Ys 


which is quite different from equation (30). 


Equations (30) and (31) were obtained by fitting the best straight lines (using the 
method of least squares) to the linear logarithmic regions close to the wall and averag- 
ing the results ahead of, and downstream of, the normal shock wave. Straight line fitting 
was not attempted when there were less than eight measured points in the linear region, 
or when the pitot traverse results contained reversed flow, The two traverse positions 


that were closest to the wall were omitted when fitting best straight lines. 


Fig 2] shows the variation of A and B with streamwise position for each 
undisturbed free stream condition, It demonstrates most clearly the reduction in A 
accompanied by an increase in B for the boundary-layer traverse positions downstream of 
the normal shock wave. There does not seem to be any consistent variation of A and B 


with Reynolds number. 


It would appear therefore that either large perturbations or normal pressure 
gradients affect the law of the wall, in a way which implies an increase in the eddy 


viscosity. 


4.5.4 Departures from equilibrium and the law of the wake 

An attempt has been made to demonstrate violent departures of the flow from equili- 
brium by extracting from the measurements the parameters used to describe the equilibrium 
locus in integral models of turbulent boundary layers. The analysis is further extended 


in terms of the character and magnitude of the wake component. 


In order to deal with separated flows, East, Smith and Mereyman’” redefined the 


Fyegey 


usual parameters 
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g = Hit]? , (32) 
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thus removing the difficulty created as the skin friction coefficient passes through zero 
and becomes negative. With an empirically derived allowance for compressibility, the 


parameters become 


— 0.04% 
Es = z (36) 
G 
— a ( + 0.04m7) : (37) 
P CG 


In Fig 22, the loci for the measured boundary layer developments are compared with 


the equilibrium locus 


Ee = 0.024 - See. (38) 


which is substantiated by experiment for attached flow and tentatively assumed in Ref 27 


to remain valid in separated flow, 


Points to the right and left of the equilibrium locus correspond respectively to 
stronger and weaker adverse pressure gradients than those appropriate to equilibrium 
flow. Negative values of ES represent favourable pressure gradients, while negative 


values of E represent separation, 


f 
In Fig 22 the results for each set of test conditions have been plotted using 


symbols to represent the following regions 


A represents equilibrium conditions ahead of the main interaction; 

B represents the strongly adverse pressure gradient under the leading 
compression (which causes separation at M = 1.5); 

C represents the rapid recovery just downstream of the normal shock wave; 

D represents the remainder of the flow which might be expected to approach 


equilibrium conditions. 


Points on the loci where there are changes from one region to the next are 


indicated by the numerals | to 3, 


ts. 


The loci show that region A is approximately in equilibrium but the rapid rise of 
. while EE only falls slowly in region B indicates that the boundary layer lags in 
its response to the strong adverse pressure gradient under the leading compression. In 
region C , Be recovers rapidly to the equilibrium locus but then overshoots before 


f 
‘his again indicates a lag in the response of the boundary layer to a change in pressure 


turning. E, only starts to increase after E has overshot the equilibrium locus and 


‘radient. The failure of the boundary layer to return to the equilibrium locus far down- 
stream (in region D) may be interpreted as a further indication of the persistently dis- 
plive effect of the shock-wave boundary~layer interactions on the velocity profiles 


viready detected in the law of the wall as shown in Fig 21. 


On the other hand the wake components of the profiles, remain fairly close to the 
«tandard shape as shown in Figs 23 and 24, In these figures the normalised wake 
v“omponents have been plotted in the following way. Fig 23a shows the undisturbed wake 
components in region A , and the results under the firsc compression through which it is 
sossible to fit a law of the wall are shown in Fig 23b. Fig 24a-c has been produced by 
arouping similar normalised wakes into one diagram. In fact Fig 24a&b covers the rapid 
recovery region C together with the beginning of D , while Fig 24c covers the down- 


stream part of region D. 


According to Winter and Gaudet? the law of the wake for compressible boundary 


layers in zero pressure gradients is 


a = 0,89]1 + sin ( = 0.483] : (39) 
UL E 0.999 
shich after normalisation and allowance for the difference between 59,999 and “995 
Ome S 
=S) = o.shi + sing ar (7 a 0.523) : (40) 
ut 4 : 0.995 


vivre the normalised wake component (av/vi) has been obtained by dividing the wake 
-omponent by the maximum value for each traverse. The shape of the normalised wake com- 
puuents has been plotted in Figs 23 and 24. On each figure a curve in the form of 


equation (40) is shown with the constants adjusted to fit the particular set of results. 


The maximum values of the wake component for each traverse (ui) are plotted 
> Fig 25 against Ee . In Fig 25a the results for regions A and B ahead of the 

sormal shock wave are given while the results downstream of the normal shock wave (in 
regions C and D) are shown in Fig 25b. It will be seen in Fig 25a that apart from 
the boundary layer profiles that have been perturbed by the leading edge of the strong 
compression, the maximum values of the wake components are in good agreement with the 
locus obtained from the equilibrium family of boundary layers described by East, Sawyer 
and mashes However, in Fig 25b for locations downstream of the disruptive effect of 

the strong pressure gradient (points in the regions C and D downstream of the normal 


shock wave) the results depart considerably from the equilibrium locus, 


“4 


ror 
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The different behaviour upstream and downstream of the shock wave is illustrated 
further in lig 26 where calculated values of J and K are plotted against the equili- 
brium paramecer EE . The values of J and K were obtained by fitting sine waves of 
the form 
le ae 0.sfi + sin (= - x} ' (41) 
ut °0.995 


aN 
to the wake components for each boundary layer profile. 


It should be noted that because of the quantity of data recorded, the values of 
(ee) were taken at that y-position of the pitot probe which gave a maximum value 
rather than estimating the true value by curve fitting. This in part accounts for the 
spread of values of J . No attempt was made to investigate the region B under the 


first compression. 


However in the undisturbed equilibrium boundary layers of region A , the averaged 


values of J and K_ result in a normalised law of the wake 


UY = 0.51 + sin 7 _ (—— - 0.494) (42) 


which is extremely close to Winter and Gaudet's law (equation (40)). The average value 
of (ens. for region A is 2,23 which is rather higher than Winter and Gaudet's 


value of 1.78 (equation (39)). 


In the region downstream of the shock wave (Fig 26b) the parameters J and K 
show an appreciable variation with Ee . Thus in this region the velocity profiles cannot 
be represented by a family having either a standard law of the wall (section 4.5.3) or a 


standard form of the law of the wake. 
5 CONCLUSIONS 


Seven flows have been studied involving the interaction of normal shock waves with 
nominally two-dimensional turbulent boundary layers over a range of Mach numbers of |.2 
to 1.5 and of Reynolds numbers based on an effective streamwise run of !10 x 10° to 
30 * 10°, The data extend from 700 mm upstream to 3000 mm downstream of the normal shock- 
wave position, These distances are equivalent to respectively 90 and 400 times the 


undisturbed boundary Layer displacement thickness. 


Oil flow investigations on the floor under the main interaction region indicate a 
highly three~dimensional flow at M = 1.5, but the results and momentum balance calcula- 
tions support the view that this three-dimensionality is confined to the separation 
region. This region is very shallow having a depth of less than !0 mm and streamwise 
extent of 300 mm while the flow in the region of the tunnel centre line is tolerably 
two-dimensional, except in the immediate regions of the saddle points at the start and 
end of the separation. Further downstream of the main interaction region, the momentum 
balance calculations indicate a slight flow divergence over the width of the tunnel of 


between 7° at M = 1.3 to 15° at M = 1.5. 


In calculating the momentum balance, use has been made of corrected static pressure 


weasurements and an attempt has been made to follow the recent ideas Ot matching tne 


Deeestream conditions with the viscous Layer by defining an equivalent inviscid flow. 


A method has been evolved for correcting the static nressure measurements made with 


‘he twin static pressure probe used in the experiment, and static pressure distributions 


rmal to the wall within the main shock-wave boundary-layer interaction are presented, 


vmbined with piitot measurements this enables the Mach number Cistribution to be vbtained. 


a : mas 3 : Senha 
fhe results agree with the laser measurements of East” which cover a more limited 
1. The present experiment has also produced more detail a: the separated regien. 


d azreement is also noted with Kooi's’ results at M = 


All the measurements confirm that the flow has not stabilised after 3 m downstream 
She normal shock wave (or 400 times the undisturbed displacement thickness}. Here the 
‘ip shear layer which is produced at the bifurcation point of the sheck waves at M = J.4 
1.5 is intact and there are still velocity gradients across the tunnel in the inviseid 
tew, Also the bou.dary-layer profiles have not recovered to the eqilibrium shapes. 


The analysis in terms of the Law of the wall and the law ct the wake for the region 
ad or the interaction agree fairly closely with Winter and Gaudet (equation (290) and 


.eS (equation (5)). However the evidence downstream of the normal shock wave sus 


one) 


the law of the wall has been changed by the large perturbations or norma! 


vadbents. However a single law of the wall may be used tor the whole regien behind the 


owtiti disturbance, 


1 


A normalised version of Winter and Gandet's law oi the wake (equation (400) provides 
‘air estimate for the shape of the normalised wake components ahead of the normal shock 
ve. The disruptive effects of the normal shock waves appear mainly as chases in 


fifude of the normalised wake cemponents, but the former the wake component fe cise 


the iirs: 


Apert trom the beundary layers immediately affected by the leading ed: 
~presslon, the maximim values ot the wake function ahead of the normal shock wave have 
imtlar .orrelation with the equilibrium function Ee (equation (36)) as the equilibrinuir 
d 
imily of flows of East, Sawyer and wash’: Although the maximum values of the wake 
saettens downstream of the normal shock wave also correlate with I. , they correlate in 


{ 
i 


quite different way from that expected for equilibrium boundary layers, 


met 


Appendix 
STATIC PRESSURE PROBE CORRECTIONS 


A.l Calibration measurements 


As stated previously (section 2.2.4) the corrections to be applied to the static 
pressure probes were based on limited measurements made in the free stream at transonic 
speeds together with the results of the calibrations of the twin static pressure prube by 
Lraverses through the floor boundary layers at free-stream Mach numbers between U.6f and 
0.88. Additional information was derived during the actual experiment from the static 
and pitot traverses ahead of the interaction region at free-stream Mach numbers of 1|.3, 
1.4 and 1.5. Data was also used from the static and pitot traverses made through the 
boundary layer 1 to 3 m downstream of the normal shock wave at M = 1.3, It was assumed 
that tor all these regions the static pressure for each traverse was constant and equal 


to the wall value, Py 


The results of the free-stream calibrations are shown in Fig 27 where the 
differences between the measured and actual static pressures are presented as pressure 
coefficients based on local conditions plotted against Mach number. Alsu (shown dotted) 
are the calibrations after three stages of smoothing which were used to correct the 
static pressure readings made under the severe velocity gradients of the region investi- 
gated (section A.2.2). There was no apparent Reynolds number effect on the calibrations 
but cime Limitations precluded much use of the single static pressure probes. fhe correc- 
tion method described, therefore applies only to the twin static pressure probe. No 


F z ; r = i) 
Static measurements were made at M = 1.4 for a unit Reynolds number of 3.5 © 10 jm, 


Awe The correction method applied to the twin static pressure probe 


Fig 28a&b shows typical static pressure measurement errors incurred while making 
boundary layer traverses with the twin static pressure probe. The boundary layers had 


constant static pressures normal to the wall. 


The errors have been presented as the differences between the measured and actual 
Static pressures and are shown as pressure coefficients based on local conditions. They 
are plotted against the distance from the wall, y , non-dimensionalised by dividing by 


the boundary layer thickness 0.995 


Fig 28a shows typical errors when the maximum Mach number is less than the critical 
value of 1,098 (the point at which the free-stream calibration errors are shown in Fig 27, 
to become negligible). The errors have a different character if the maximum Mach number 
tor the traverse is above the critical value as shown in Fig 28b. It is therefore con- 


venient to divide the errors to be corrected into two categories: 


(1) when the outer static pressure probe is at a Mach number of less than 


or equal to 1.098; 


(2) when the outer static pressure probe is at a Mach number greater than 


1.098, 


This is a not unreasonable assumption because the main source ot error in the tree 
ream calibrations and the statie traverse measurements made away trom the int lies. 
tie wall, can be shown to be mostly due to the interconnecting stem between the inner and 
) 

“ter static pressure probes (see Fig 27 and Revd. Pope and Cooksey’, pp loséloge. 
serefore the main correction to be applied to the inner static probe depends on som 
estimate of the mean Mach number over the connecting stem provided the static pressure 
“robes are both at Mach numbers below the critical |.098. If however, the outer static 
‘rebe is at a Mach number above this critical value, then there is no interference 
-cseciated with the interconnecting stem in the region of the outer static probe. In 

is case a suitable course of action is to use the local Mach number to estimate the 


.<rrection due to the interconnecting stem on the pressure measured by the inner static 


vrebe, 


it is now proposed to deal with the errors in more detail by considering the three 
vesions A, B and C indicated in Fig 28 in conjunction with the Mach number of the 
eler static pressure probe. 


\.2.t Region A 
The error in this region is treated as arising mainly from wall interference. There 
However some compressibility effect. The errors are presented in law of the wal! 
vrms for correlation purposes (Figs 29 and 30). Fig 29 shows the accumulated errors 
resulting from the subsonic boundary-layer traverses, while Fig 30 shows the errors 
“csulting from the undisturbed boundary-layer traverses ahead of the interaction region at 
tree-stream Mach numbers of 1.3, 1.4 and 1.5. To calculate the corrections, the errors 
region A may be added to those of region B . However because the errers of regions 


5 and C are based on local conditions rather than the conditions at the edge of the 


“sandary layer used for region A , the region A errors are multiplied by a factor 


ae S| to bring them into line. 
Thus the corrections needed for the static pressure measurements made in region A 


given by 
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The measurement errors shown in Fig 29 have been pres 


Pa Pw) 
2qyc 
do Uy f 
and plotted against yl /¥, . The corrections needed are actually fairly small ‘fer the 


Mach numbers are low) and so a single correction has been used and this is represented 


by the continuous line in the figure. The correction is given by 
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= isc] when ag < 644 
A = 
5) 
1.\7 UL yu, 
C = c.l— 50.37 - 5.469 Inf— when 644 <— < 1000 
Pp f£\(M ve ae 
A § 
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c = 0 when —— 2 1000 
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A.2.1.2 Outer static pressure probe at M > 1.098 


Fig 30 shows the pressure measurement errors plotted in similar terms to those 
used in section A.2,1.1. Here the results are strongly Mach-number dependent. The 
corrections have been applied on the assumption that the pressure measurement errors in 
Fig 30 have a linear relation with Mach number for constant yl_/., . The slopes and 
zero intercepts of this linear variation are then assumed to be entirely dependent upon 


yU_/v¢ The corrections which are shown by the continuous lines in Fig 30 are 


0.995 ° 
given by 
M, : ye 
Cc of re) M = 0.2725 - 0.1153 in —) 
Pa I\M Nt 
yu yu. 
where Fo = 19,33 inf) - 433.5 when —— < 7660 
8 "ss 
yu yl, 
F = 692.7 inf ) - 6456 when 7660 <—— < 11200 
& eR 
yU_ 
F = 0 when —— 2 11200 


This correction is based on rather limited evidence but the magnitude is small 


and the probable accuracy of the correction is within 14% of the true static pressure. 


A.2.2 Regions B and C 


The errors in these regions are treated as arising from a combination of shear 
flow and compressibility effects. It would therefore be expected that there would be an 
error which would correlate with y/* 9 995 and a larger error due to the geometry of the 
twin static pressure probe, which would correlate with Mach number. As stated in section 
A.2 the greater part of the latter error is due to the stem interconnecting the inner and 
outer static pressure probes, and the necessary correction depends on an estimate of Mach 


number and a form of the free-stream calibration of the twin static probe. 


In fact, the average Mach number of the inner and outer static probes is used te 
give the necessary correction to the static pressure measured by the inner probe, when 
the outer probe is below the critical Mach number of 1.098, In all other instances the 


local Mach number is used to obtain the correction needed to the measured static pressure. 
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This applies to both the inner and outer static probes. The local Mach number is always 
used for the outer static pressure probe because it is normally in a region where the Mach 


number gradient normal to the wail is small. 


It was found that the free-stream calibration overcorrected the static pressure 
Measurement errors, and so a compromise calibration was used (shown dotted in Fig 27a) 
which involved three stages of smoothing. It may be noted from Fig l6a~-f that the criti-~ 
cal Mach number is only encountered ahead of and in the immediate vicinity of the main 
interaction region. The location of the point for critical Mach number is close to the 
wall for stations upstream of the influence of the normal shock wave. The environment is 
therefore one of severe shear strains and velocity gradients and the smoothing exercise 
is probably justified as it only affects the free-stream calibration in the region of the 


eritical Mach number. 


The form of the corrections is given by 


Cc en(y x) ; 


0.995 


A.2.2.1 Outer static pressure probe at M< 1.098 


The corrections are calculated in terms of local density and velocity. The modi- 
2 
tied calibration curves (shown dotted in Fig 27) are used to obtain (Pp 7 P/ acu") from 
fn(M) and the modification to the free-stream calibration due to probe position is shown 


in Fig 31. 


Thus the total correction is given by 


M + ? , f 
Cc = fn (2) x 0.94] 8in{0.665 (2 - 0.15] 
Px “0.995 \ 


fon 


C = fn(M,) * 0.94 vino. (<- - 0.15) 
Pc 0.995 

c = 0 if —Y- < o.15 . 
Pa 0.995 


4.2.2.2 Outer static pressure probe at M > 1.098 


& 


c = fn(M,) » 0.94 rin}0.665 (<— - 0.15) 
0.995 


x 


Appendix 
} 
C= fn (,) + 0.94 vin fo.665(—t— - 0.13)| 
Pc 0.995 
G0 if —— < u.15 . 
PR 0,695 
A.2.3 Total corrections 
ee 
Inner tube C = C #5€ = > 
I Pa Pa keen 
Po- p 
Outer tube Cc eo = m = 7 
PL Pc beue 


A.3 Applying the corrections 
The following describes how the static pressure corrections are calculated and 
applied: 
(1) Calculate the Mach number at each measuring station from interpolated pitot and 
wall pressure measurements. Calculate ylig 995° Ce and y(U_/v,). 
C2) Calculate the static pressure corrections and subtract from the measured values. 
3) Recalculate the new Mach number using the new static pressure and the interpolated 
pitot measurement. 
Steps (2) and (3) are repeated until the difference between the new and old Mach 
number is less than 0.001. A suitable damping factor is needed to prevent oscillations 
in the calculation, 


(4) 


Record the corrected Mach number and static pressure. 


Fig 32 shows the effect of the correction on static pressures measured ahead and 


downstream of the interaction region where the static pressures may be expected to be 


nearly 


uniform across the boundary layer thickness. It will be seen that the static 


pressure errors have been reduced to 14% of the true static pressure. 


Aw4 


Transferring the corrected static pressures to the pitot traverse measurements 
a eg a EE OE NEE SE MC OSUL SMES 


The final stage of the procedure is to interpolate in the distribution of static 


pressure to obtain values at the location of the pitot-pressure measurements. 


The calculated static pressures for each static pressure traverse are smoothed and 


transferred to a carpet plot of static pressure against streamwise position as show in 


Fig 33. 


The appropriate static pressures are then transferred from the carpet plot for 


each pitot traverse, A final stage of smoothing ensures that the final static pressure 


curves 


floor. 


are asymptotic to the wall static pressure within a few millimetres of the tunnel 


The corrected static pressure results may be seen in Fig !5a&b. 
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LIST OF SYMBOLS 


(Bracketed symbols represent computer headings) 


‘law of the wall' constants in equation (28) 


local skin friction coefficient obtained from East's prediction 
of tha law of the wall; see equation (6) 


local skin friction coefficient based on Patel's Preston tube 
calibration as formulated by Head and Vasanta Ram; see section 3.3 


pressure coefficient based on local conditions 

probe diameter 

constant in equation (6) 

equilibrium factor in equation (36) 

equilibrium factor in equation (37) 

compressibility factor in equation (6) 

Clauser shape parameter in equation (32) 

distance of probe centre line above wall 

boundary layer shape parameter in equations (9) and (13) 


boundary layer s*.1: + parameter in equations (10), (14) and (20) 


6 - 6% 
boundary layer shape parameter 0-385 


normalised ‘law of the wake’ constants in equation (41) 
effective turbulent run ahead of normal part of shock wave 
Mach number 

Mach number at edge of boundary layer 

static pressure 

equivalent inviscid static pressure 

tunnel total pressure 

local total pressure 

unit Reynolds number; see equation (3) 

temperature in Kelvins 

total temperature in Kelvins 

horizontal component of velocity 

horizontal component of velocity at edge of boundary layer 
wake component of velocity profile 


friction velocity 
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LIST OF SYMBOLS (concluded) 


ut equivalent incompressible friction velocity; see equation (26) 
T 
Vv vertical component of velocity 
x streamwise position relative co the normal part of the main siock 
wave 
y CY) equivalent height above wall; see equations (1) and (2) 
y* non-dimensionalised y with respect to wall parameters in 


equation (6) 


Zz distance perpendicular to flow but parallel to surface 
& boundary layer thickness 
59 995 (DEL995) boundary layer thickness at U/U, = 0.995 
é boundary layer thickness at U/U, = 0.999 
0.999 y aay /U, 
‘e (DELE ) boundary layer energy thickness; see equation (15) 
5 (DELBAR) bounda.y layer thickness parameters; see equation (7) 
hal (DEL*) boundary layer displacement thickness; see equations (11) and -!%: 
§ (THETA) boundary layer momentum thickness; see equations (12) and 4 
a (THETABAR) boundary layer thickness parameter; see equation (8) 
K mixing length constant in equation (27) 
“ surface curvature 
viscosity ot fluid 
kinematic viscosity of fluid 
i} pressure gradient parameter for equilibrium flows; see equaticn «+ 
Gensity 
Subseripts 
i eguivalent inviseid tlow quantities 
1 inner static pressure probe 
m measured quatities 
normalised quantities 
0 tree-stream conditions betore start of interaction 
Pp intermediate equivalent inviscid tlow quantities 
U outer Static pressure probe 
w wall conditions 
A,B,C regians in Fig 28a&b 
‘ edge of boundary layer conditions 
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Fig6b Schlieren photographs of normal shockwave boundary-layer interaction 
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Fig9 Momentum balance based on measured p + calculated p; for X < 0 
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Fig 10a Flow divergence needed to balance momentum integral equation, M = 1.3 
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Fig 10b Flow divergence needed to baiance momentum integral equation, M = 1.4 
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Comparison of momentum balance calculations at Re/m = 10 x 10° 
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Fig 13b Boundary layer development (p = Py) M = 1.4 
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Fig 13c (concld) Boundary layer development (p = Py) M = 1.5 
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Fig 18 — Logarithmic velocity profiles ahead of interaction (region A) (see Fig 22) 
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Fig 26a%b Constants J and K in (AU/U!)qy = 0.5(1 + sin 7/Jly/59 995 — K)) 
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Fig 27a&b Free stream calibration of static probes 
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a) Behind shockwave 
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b) Ahead of shockwave 
X=-523mm Mp5 = 1.284 


Fig 28a&b Typical static pressure measurement errors 
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Fig30 —Wali/static interference errors, My 2 1.098 
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Fig 31 
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Fig 31 Modification to smoothed free stream static pressure calibration caused 
by probe rosition 
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Fig 32a&b 
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Fig 32a&b Typical corrected static pressures 


Fig 33 


TA 92099 


Fig 33 Smoothed static pressure distribution, M = 1.5, Re = 10x 106/mm 
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